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Abstract 

Massive MIMO and small cell are both recognized as the key technologies for the future 5G wireless systems. 
In this paper, we investigate the problem of user association in a heterogeneous network (HetNet) with massive 
MIMO and small cells, where the macro base station (BS) is equipped with a massive MIMO and the picocell BS’s 
are equipped with regular MIMOs. We first develop centralized user association algorithms with proven optimality, 
considering various objectives such as rate maximization, proportional fairness, and joint user association and resource 
allocation. We then model the massive MIMO HetNet as a repeated game, which leads to distributed user association 
algorithms with proven convergence to the Nash Equilibrium (NE). We demonstrate the efficacy of these optimal 
schemes by comparison with several greedy algorithms through simulations. 

Index Terms 

Massive MIMO; small cells; heterogeneous networks (HetNet); user association; unimodularity; game theory. 


I. Introduction 

Over the past two decades. Multiple Input Multiple Output (MIMO) has evolved from a pure theory to a practical 
technology, and has greatly enhanced the wireless system capacity by offering many degrees of freedom (DoF) 
for wireless transmissions. However, due to the so-called “smartphone” revolution, mobile users are demanding 
increasingly higher data rates for rich multimedia applications. Existing and future wireless networks are facing the 
grand challenge of a 1000-time increase in mobile data in the near future JT). There have been tremendous efforts 
made aiming to cater for this demand. For example, based on MIMO and OFDM, LTE-Advanced targets at a peak 
rate of 1 Gbps, but the average rate is still less than 100 Mbps. In the foreseeable future, such rates can hardly be 
satisfactory for data-hungry wireless users. 

To boost wireless capacity, two technologies have gained most attention from both industry and academia. The 
first one is massive MIMO (a.k.a., large-scale MIMO, full-dimension MIMO, or hyper MIMO) J2], 0. The idea 
is to equip a base station (BS) with hundreds, thousands, or even tens of thousands of antennas, hereby providing 
an unprecedented level of DoF for mobile users. The massive MIMO concept has been successfully demonstrated 
in recent works a, 0 . The second technology is small cell. A great benefit of deploying small cells is that 
the distance of the user-BS link can be effectively reduced, leading to reduced transmit power, higher data rate, 
enhanced coverage, and better spatial reuse of spectrum. Both massive MIMO and small cells are recognized as 
key technologies of the future 5G wireless systems g). 

In this paper, we consider a heterogeneous network (HetNet) with massive MIMO and small cells, where the 
macrocell BS (MBS) is equipped with a massive MIMO and the picocell BS’s (PBS) are equipped with regular 
MIMOs. To fully harvest the benefits promised by these two technologies in an integrated HetNet system, it is critical 
to investigate the user association problem, i.e., how to assign active users to the BS’s such that the system-wide 
capacity can be maximized and users’ experience can be enhanced. 

There are already several recent works pushing forward in this direction. In IITI- flOl . the authors consider the 
problem of user association in massive MIMO systems operated in the frequency-division duplexing (FDD) mode. 
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These papers are focused on a macrocell without small cells. In fTTI . user association in time-division duplexing 
(TDD) massive MIMO system is addressed, where factional user association is allowed. Bayat et al. in fl2l model 
the problem of user association in a femtocell HetNet as a dynamic matching game and derive the optimal user 
association. However, massive MIMO is not considered in the system model. In fl3l . the authors investigate the 
problem of user association with conventional MIMO BS’s and propose a simple bias based selection criterion to 
approximate more complex selection rules. Bjornson, et al. in lfl4l consider the problem of improving the energy 
efficiency without sacrificing the quality of service (QoS) of users in a massive MIMO and small cell HetNet. 

Motivated by these interesting works, we consider the user association problem in a TDD massive MIMO HetNet 
in this paper, taking into consideration of the practical constraints, such as the limited load capacity at each BS, while 
without allowing fractional user association. The main goal is to maximize the system capacity while enhancing 
user experience. 

More specifically, this paper contains two parts: (i) centralized user association and (ii) distributed user as¬ 
sociation. For centralized user association, we investigate the problems of rate maximization, rate maximization 
with proportional fairness, and joint resource allocation and user association. We prove the unimodularity of our 
formulated problem and develop optimal user association algorithms to the problems of rate maximization and rate 
maximization with proportional fairness. We then propose a series of primal decomposition and dual decomposition 
algorithms to solve the problem of joint resource allocation and user association and prove the optimality of the 
proposed scheme. For distributed user association, we model the behavior and interaction between the service 
provider, who owns the BS’s, and users as repeated games. We consider two types of operations: (i) the service 
provider sets the price and the users decide which BS to connect to, and (ii) the users bid for the opportunity of 
connection. We prove that in both cases the the proposed algorithms converge to the respective Nash Equilibrium 
(NE). 

In the reminder of this paper. Section |TT] introduces the system model and preliminaries. Optimal centralized 
and distributed user association schemes are presented in Sections [III] and IIVI respectively. Section [V] presents the 
simulation study and Section [VI] concludes this paper. Throughout this paper, we use a boldface upper (lower) case 
symbol to denote a matrix (vector), and a normal symbol to denote a scalar. (•) H denotes the Hermitian of a matrix. 


II. System Model and Preliminaries 

The system considered in this paper includes K users and J BS’s, including an MBS with a massive MIMO 
and (J — 1) PBS’s, each equipped with a conventional MIMO. The channel model is hj^, n = gj,k,nlj,k, where 
hj t k,n is the channel of antenna n at BS j to user k, fjj.k.r, represents the small scale fading coefficient between 
antenna n of BS j and user k, and l :l j ; stands for the large scale fading coefficient between BS j and user k d. 
Concatenating all the channel coefficients from all the antennas of BS j , we obtain the channel vector hj.fc, as well 
as the channel coefficient matrix for signals transmitted from BS j as H, = [hj i, hj 2 > ■ • ■ , liyj, : ] . 

Let y-j denote the signals received by the users connecting to BS j, W, the precoding matrix of BS j, and d y 
the data sent from BS j. We have 

y, - H ; w,(l, - n ( . (1) 

where n j is the zero mean circulant symmetric complex Gaussian noise vector. 

Each active user has the options to connect to either the MBS or a PBS. For a user k, define user association 
index variable Xk :i as 

{ 1, if user k is connected to BS j. 

( 2 ) 

0, otherwise. 

Let its achievable rate if connected to BS j be Rk r r)^ = XkjRk r and its actual data rate be rjk- We have 

( 3 ) 

3 3 

For users connecting to a massive MIMO BS j (i.e., the MBS), their achievable rate can be approximated with 
the following deterministic rate HD- 


„ , / M, -Lj + 1 

Rk, = log l + ——- 1 - 
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j l j,k 




( 4 ) 



where Mj is the number of antennas at the BS, Lj is the prefixed load parameter of the BS indicating how many 
users it could serve, and Pj is transmit power from the MBS. Note that there is no small scale fading factor in (f4]>. 
This approximation has been proven to be accurate hd. 

For a PBS with a conventional MIMO, we assume that the inter-cell interference is negligible among the picocells, 
due to the small transmission powers and effective inter-cell interference coordination (ICIC) 1131 . The achievable 
rate of user k connecting to PBS j can be represented as follows. 


Rkj = log 
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(5) 


where w j tk is the fc-th column of BS j’s precoding matrix W,. There are many precoding designs for conventional 
MIMO BS’s, such as matched filter (MF) precoding, zero forcing (ZF) precoding, and regularized zero forcing 
(RZF) precoding (16). Without loss of generality, we adopt MF precoding in this paper with W j = -l===Flj / , where 
ip is a power normalization factor. The signal received by all the users connecting to PBS j can be rewritten as 
follows. 


^ h^hjidi + h^ 1 hj^d2 + • • • + hj-jij^dk^ 
hfphjpdi + hj^dz + • • • + hR^hjykdk 


( 6 ) 


\ h f,k h i, + h f,k h j,2d-2 + • • ■ + hf tk h jt kd k 


Thus, the achievable rate for user k regarding to PBS j can be obtained as follows. 


Vkj = log 
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III. Centralized User Association 


In this section, we consider the problem of centralized user association. We assume that the BS’s have all the 
channel state information (CSI) via uplink training. We adopt the following utility function for each user k with 
achievable rate r] k . 


U(jq k ) = < 


ri\~ a /(l-a)i ifa> 0 ,a^l 
Vk, if a = 0 

log ( 77 k), if ol = 1 . 


( 8 ) 


When a = 0, maximizing U(-) yields the maximization of the sum rate (but no fairness); when a —> oo, it leads 
to the maximization of the worst-case rate (i.e, max-min fairness); when a = 1, it yields the maximization of the 
geometric mean rate (i.e., proportional fairness). 

Our goal is to maximize the system utility by configuring the user-BS association. Typically, we consider the 
cases when a = 0 and a = 1. In the case of a = 1, we define U( 0) = 0. 


















A. Maximizing Sum-rate 

We firstly investigate the problem of maximizing the system sum rate, i.e., a = 0 in © and U{r]k) = r) k . The 
problem can be formulated as follows. 

K 

Pl-1: max rjk (9) 

s.t. 'Yhx kj < Lj < Mj, j = 1,2, • • • , J 

k 

y, x k j <1, fc = 1,2,--- , K 

3 

Constraints ©, 0, (Q}. ©. 


Note that the second constraint requires the number of users connecting to a BS to be no more than its prefixed 
load, which should in turn be no more than the number of antennas it has, since theoretically BS j can provide 
at most Mj degrees of freedom (DoF). Assuming the Lj’ s are already chosen to satisfy Lj < Mj, we drop this 
constraint in the remainder of this paper. The third constraint simply claims that each user can connect to at most 
one BS at a time. 

A key observation is that ([7} can be rewritten as 


V kj = x kj log 

Thus the R kj in ([U can be redefined as 

R kj = log 
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( 11 ) 


In (im >, it can be seen that R kj depends on other users’ choices x k >., for all k ^ as well. To make the problem 
tractable, we adopt the worst-case approximation by assuming the users within the coverage of BS j (denoted as 
Qj) all connect to BS j with perfect channels. This way, O can be approximated as 


R kj = log 


1 + 


Pj 


2 

i + (\gj\-i)Pj 


( 12 ) 


where j | for a set stands for the cardinality of the set. 
Define auxiliary variables c kj as follows. 


° k j = 


R kj in ©, if BS j is the MBS; 
R kj in O, if BS j is a PBS. 


(13) 


















The sum rate maximization problem can be reformulated as 


K j 

Pl-2: max EE XkjCkj (14) 

' Xk i i fc=l j=1 

s.t. 'y ' Xkj ^ Lj, j = 1,2, • • • , J 

k 

E 2 ^ < k = 1,2, ■ ■ ■ ,K 
j 

Constraints ©. (fl3l ). 

Since the variables Xk, ’s are binary, problem Pl-2 falls into the category of Multiple Knapsack Problems, which 
is one of Karp’s 21 NP-complete problems uni- Although a greedy algorithm could be developed to compute 
sub-optimal solutions, we show that problem Pl-2 can actually be optimally solved by taking advantage of its 
special structure. 

Let X be a matrix with entries Xk r k = 1, 2, ■ • ■ ,K,j = 1, 2, ■ • • , J. We could convert X to a vector x by 
concatenating the rows of X and taking a transpose as x = [xi x X 2 , ■ ■ ■ Xk, ■ ■ • X\ s ■ ■ ■ Xkj] T , and simplify 
the notation as x = [xi X 2 • • • xkj] ■ We then apply the same conversion to the matrix comprising Cfc. and obtain 
vector c. Problem Pl-2 can be rewritten as 

Pl-3: maxc T x (15) 

X 

K 

s.t. y ^ ‘ Lj, j = 1, 2, • ■ ■ , J 

k=l 

J 

y \ Xk~\ -(j—i)K — k 1,2, • * • , K 
1=1 

Constraints Q. ( |T3T i. 

Ignoring constraints © and (fl3l >. dehne A as the constraint matrix of problem Pl-3, with entries being the 
coefficients of the first and second constraints. We next introduce an important definition and derive a key lemma. 

Definition 1. A matrix A is called totally unimodular if the determinant of every square submatrix of A is either 
0, +1 or -1 /[7g| J . 

Lemma 1. The constraint matrix A of problem Pl-3 is totally unimodular. 

Proof: Inspecting the constraints in problem Pl-3, we find that A is of the following form. 
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( 16 ) 




We can divide A into blocks as follows. 


A 


Ai A 2 • • • A j 
Bi B 2 ■ ■ ■ B j 


(17) 


where each A j, j £ [1, J], is a submatrix of A of size J x K; and each B y . j £ [1, J], is an identity matrix of 
size K x K. 

Let S n denote an arbitrary square submatrix of matrix A of size n. For any submatrix of A of size n = 1, it is 
trivial to see that the determinant of this submatrix is either 0 or +1. So we only need to consider the case where 
the size of the square submatrix is greater than or equal to 2. 

Case 1: S n is taken entirely from one of the submatrices A, or B y , j £ [1, J], We can see from the structure 
that at least one row of A j is all zero. So if the square submatrix is entirely taken from A,, the determinant of the 
submatrix is zero. Since matrix B j, for all j, is simply an identity matrix, it is straightforward that the determinant 
of any square submatrix of B j is either 0 or +1. 

Case 2: S„ is not entirely taken from any one of the submatrices A j or B ( , j £ [1, J], In this case, the square 

submatrix must be taken from 2n (n = 1, • • • , J) submatrices of the submatrix set (A j U Bj, j £ 1, • • • , J). We 

next proceed with our proof by applying induction method. 

For the base case n = 1, the square submatrix to be examined is of size 2. Since the entries can only be 0 or 
+1, the determinant can only be 0, +1 or —1. 

Now assuming that any square submatrix of size (n — 1) has determinant 0, +1 or —1, we need to check if the 
same conclusion holds for any square submatrix of size n. 

We first notice that each column of A has exactly two +ls. Moreover, exactly one of them is in A j, and the 
other in Bj. Let q* = argmin ? JL S ni q , where S ni q is the (i,q)- th entry of S„. That is, column q* has the 
minimum number of Is among all the columns of S„. 

Let Qq* = min g S rai . can only be 0, 1, or 2. 

If C/* = 0, then all the entries of the r/'-th column of S n are 0, which results in det( S„) = 0, where (let, is short 

for determinant. 

If (q. = 1, then we could calculate det(S ra ) by expanding the q*-th column and obtain det(S„) = det(S( ra _ 1 )). 
Since det(S( n _ 1 )) is 0, 1 or —1 by our induction hypothesis, we conclude det(S„) is 0, 1 or —1. 

If = 2, we could firstly negate all the entries taken from Bj, and then add all the rows in Bj to any non-zero 
row in A } . After this procedure, if that non-zero row in A j is still non-zero, add that row to any other non-zero 
row in A j. Repeat this process until we get a zero row in A j. The reason why this process always give us a 
all-zero row is that we have equal number of +ls in A j and Bj. Since any basic row operation does not change 
the determinant and we finally get a all-zero row, we have det(S„) = 0. That completes our induction. ■ 

Fact 1. For a linear programming problem, if its constraint matrix satisfies totally unimodularity, then its has all 
integral vertex solutions H / 

Fact 2. For a linear programming problem, if it has feasible optimal solutions, then at least one of them occurs at 
a vertex of the polyhedron define by its constraints m 

Given the facts and Lemma [Q we have the following theorem. The proof is straightforward and omitted. 

Theorem 1. The optimal solution of problem PI can be obtained by solving a relaxed problem where the variables 
Xkj are allowed to take real values between [0,1], 

Given the above theorem, we could obtain the optimal solution of PI by solving the relaxed problem, termed 
NP1, using common LP solvers llT8ll . 

B. Proportional Fairness 

In this section, we take proportional fairness among user achievable rates into consideration. The problem can 
be formulated as follows. 



( 18 ) 


s.t. same constraints as problem Pl-2. 


Problem P2-1 is a nonlinear integer programming problem, which is generally NP-hard. To get a better under¬ 
standing of the problem, we examine its equivalent problem as follows. 
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P2_2: “ ax II E 




(19) 


3 fc=i \j=i 

s.t. same constraints as problem Pl-2. 


Problem P2-2 is a geometric programming problem, with binary variables. The objective function is a posynomial 
function with J K terms. Conventionally, to solve geometric programming problems we need to introduce new 
variables such as y = log(tc) so that geometric programming can be solved via convex programming. However, 
here x kj ’s are binary. Since log(O) = -oc, we could not apply these techniques. Another heuristic scheme is to 
firstly sort these J K coefficients, and then find L :i maximal coefficients for each BS. However, even sorting these 
J K coefficients could be computationally prohibitive even for a small system, which requires 0{ J K log( J K )) 
operations. 

A key observation about the logarithm function is that log(E; T i) < Ei log(r<), for all Ti ^ 2. Therefore, in 
practice^ the optimal value of problem P2-1 is upper bounded by that of the following problem. 


K j 


NP2: max EE x kj log(Cfc 3 .) 


( 20 ) 


3 fc=l j =l 

s.t. same constraints as problem Pl-2. 

We have the following results for the transformed problems. 

Lemma 2. Problems P2-1 and NP2 are equivalent. 

Proof: Recall that if t] k = 0, we define U(rj k ) = 0. The second constraint E'E x k, < 1 imposes that 
each user could only connect to one BS. Consequently, E,-%l°g( c b) = l°g(E i x fe c fe )■ Furthermore, we have 

Efe Ej x k 3 log (c kj ) = E k lo g(Ej x k 3 c kj ). ■ 

Comparing problems NP2 to Pl-2, we find they are actually equivalent. Thus we can obtain the optimal value 
of P2-1 by applying the same technique used to solve problem Pl-2. We then have the following lemma. 


Lemma 3. Sum rate maximization in Section \III-A\ also achieves proportional fairness. 

Lemma 4. The optimal value of problem P2-1 is upper bounded by UB\ = Ea-Li log 2 (cfc J ). 
Proof: Denote m = max{ln(cfc 1 ), ln(cfc 2 ), • ■ • , ln)^)}, we have 


log 2 



< log 2 E 




ln(cfc,-) 

e m 


= log 2 (e m ) + log 2 

< mlog 2 (e)+log 2 (J) (21) 

The hrst inequality is because x k < 1. The second inequality is due to the fact that rri is the largest one among 

all the x kj c kj and e ln ( Xk j Ck f- rn < \ 

On the other hand, it follows the constraint Ej_i x k 3 < 1 that log 2 ^Ej = i x k 3 c kj ^ < log 2 (e m ), which is a 
better bound than (OH . We thus have UBi = Ea~i max j l°g2( c fcj)- ■ 


'Recall that c^. is the achievable rate of user k connecting to BS j. c^. > 2 is generally satisfied in current wireless systems with a 
sufficiently large bandwidth and high transmission power. 




Algorithm 1: Greedy Algorithm 1 for User Association 


1 Initialize K, = {1, 2, • • • , A'}, Lj,Vj € J and Xk :j to be an all-zero matrix ; 

2 for k = 1 to K do 

3 for j — 1 to J do 

4 | Compute Ckj as in O ; 

5 end 

6 end 

7 while 3j, Lj / 0 do 


8 
9 
10 
11 
12 

13 

14 end 


Find ( k*,j *) = argmax fcj {cfcj} ; 

if Lj* 0 then 

= l ; 

3 

Lj* = Lj* — 1 ; 

K. = K\k * ; 

end 


Algorithm 2: Greedy Algorithm 2 for User Association 


1 Initialize A2 = {1,2, 

2 for k = 1 to K do 


, A}, Lj, Vj £ J7 and to be an all-zero matrix ; 


3 for j = 1 to J do 

4 | Compute Ckj as in 03 : 

5 end 

6 end 

7 for j = 1 to J do 


8 
9 
10 
11 
12 

13 

14 end 


while Lj ^ 0 do 

Find (k*,j) = argmax fe {c fc ,j} ; 
x k * = 1 ; 

3 

Lj = Lj 1 ; 

K, = K\k* ; 

end 


For comparison purpose, we propose two sub-optimal greedy algorithms, i.e.. Algorithms |T] and [3 as bench¬ 
marks. They can be directly used for comparison with problem Pl-1. To compare with problem P2-1, in Al¬ 
gorithm [Q we need to change Steps 7 and 8 as “while 3 j, Lj ^ 0 & max^, log(cfc,j) > 0 do” and “Find 
( k*,j *) = argmax fe ^{log(cfc ? j)},” respectively. In Algorithm^ we need to change Step 8 and 9 as “while Lj ^ 0 
& max fe log(cfcj) > 0” and “Find ( k*,j) = argmax fc log(cfcj),” respectively. 














C. Joint Resource Allocation and User Association 


In this section, we take resource allocation into account. Consider a massive MIMO OFDMA HetNet. In OFDMA 
systems, such as LTE, the time-frequency resource is divided into resource blocks (RB). A typical RB consists of 
12 subcarriers (180kHz) in the frequency domain and 7 OFDMA symbols in the time domain (0.5 ms). So the 
system may have up to several hundreds of RBs. We normalize it to be a unit number. A user k connecting to a BS 
j gets a portion /3 kj of the overall resource. The goal is to maximize the system utility considering both resource 
allocation and user association. 

Considering the logarithm rate utility and defining <f>j = {k \ x k . = 1}, the problem is formulated as follows. 

* l 

P3-1: max ) log 

v 

s.t. ^ < 1, j = l,2 ,--- ,J 
same constraints as problem Pl-2. 

To solve problem P3-1, we need to: (i) select users for each BS to serve and (ii) allocate resources to the associated 
users at each BS. We next propose a series of primal decomposition and dual decomposition to solve the problem 
optimally. 

It is worth noting that the problem can also be formulated in a different way, by substituting constraint J2j=i x k, < 

1, k = 1, 2, • • • . K with a new constraint x k s = 1, k = 1,2, • • • , K. We call this problem P3-2. Comparing 
these two formulations, we have the following observations. 

1) Problem P3-1 does not require that every user must be connected, while problem P3-2 requires each user be 
connected, even under unfavorable conditions. 

2) Problem P3-2 has a more stringent requirement than problem P3-1. Therefore the optimal value of problem 
P3-2 is upper bounded by that of problem P3-1. 

3) Since problem P3-1 offers more choices of user association, problem P3-1 is slower in convergence than 
problem P3-2. 

We focus on the harder problem P3-1. Given the algorithm to solve problem P3-1, problem P3-2 can be readily 
solved. Due to integer variables x kj and real variables fi k . , problem P3-1 is a mixed integer nonlinear programming 
problem (MINLP), which is generally NP-hard. However, next we propose an algorithm to obtain its optimal solution. 

Since x kj ’s take binary values and J2j=i x k 3 <h we have J2k =l lo S (Z)/=i x k } c kj Pk 3 ) = EfcLi 2j=i x k 3 log(% Pk 3 )■ 

Recall that if ]Cj=i x k 3 = 0- the logarithmic utility is 0. Thus problem P3-1 can be reformulated as 

K j 

P3-3: max Y' Y' x k Tog(c fc j3 k ) (23) 

{x k .,B k .} t-—' 

\ *, r fc =1 j =1 

s.t. same constraints as problem P3-1. 

The choices of fj kj rely on the values of x kj . Given these coupled variables, we first apply the Primal Decom¬ 
position method 1201 to decompose problem P3-3 to the following two levels of problems. Fixing variables x kj ’s, 
we have the lower level problem as 

K j 

Xk :i l0 SK ) (24) 

s.t. 5Z Pkj < = M,--- ,j. 


j =1 


c ki Pk, 


( 22 ) 


When the /3fc/s are fixed, the higher level problem (or, the master problem) is given by 

K j 

max EE Xkj \og{c kj /3 kj ) (25) 

{xk > } k=ij=i 

s.t. same constraints as problem Pl-2. 


Since there are no couplings among the subproblems, the lower level problem (l24l) can be further decomposed into 
L subproblems as follows. 


K 


max E^ fe i lo s( Cfe A') 


k =l 


(26) 


K 


S.t. 


E fa, < i,i = i, 2 , - - - ,j. 


ke®i 


Defining Lagrange multiplier A, the Lagrangian of problem ( f26b is defined 


K 


K 


c = E x kMk*ifci) + A 1 - E fci ■ 


*:=i 


fc=i 


Applying KKT conditions |[2TI . the optimal solution can be obtained as follows. 

Xk; 


Substituting 


Z^fc=l X kj 

into the master problem, the objective function becomes 

K j 

EE^ 10 ® 




, Z_/fc=l X kj 


(27) 


(28) 


(29) 


fc=i j=i 

Note that we have dropped one term in ( 1291 . since due to the definition Q, we have )' 2 = Xk r Since 

Uk=i Xkj is in the denominator, problem ( |29[ ) has coupled objectives. The main idea of addressing the coupled 
objective is to introduce auxiliary variables and additional equality constraints so that the coupling in the objective 
function is transferred to coupling in the constraint l20l . We thus introduce a new variable, which is defined as: 


K 

= ^2, Xk i ■ 

fc= 1 


(30) 


To solve the above problem, we relax Xkj to a real number in [0,1], However, we will show later that even if we 
have relaxed the variables, we could still find the optimal solution to the original problem. The relaxed problem to 
be solved is 


max 


E E “fcjks (% 

fc=i i=l 


(31) 


s.t. H j < Lj, j = 1,2, • • • , J 

j 

E Xkj <1, k = 1,2, , K 

j=i 

0 < Xkj < l,for all k,j 
Constraints (IT3l) . (f30t . 


Problem ( |3TI ) is a convex optimization problem. Defining Lagrange multipliers for the equality constraints ( l30l ). 
problem (fTil) can be solved with the dual decomposition method. Alternatively, we propose Algorithm [3] to obtain 




the optimal solution of problem ( |3H j9l , 1 1 1 01 , l22l . In Algorithm^ 6® is the step size at the i-th iteration given 
by 


where •& and 7 are positive numbers. 


6® 


t + 7 ’ 


(32) 


Theorem 2. Algorithm [7] optimally solves problem ( U7I ). 


Proof: Let x® denote the solution produced by Algorithm [3] at step t. Let dU(x ®) be the subgradient of 
the objective function in problem (OH at step t. It can be easily verified that the updated direction in step [16] of 
Algorithm [3] is the subgradient direction. Since Sj is upper bounded by Lj and I \, and Xk :i is upper bounded 

by I\, dU(x®) ' s a l so bounded. 

Denote U a as the final result produced by Algorithm [3] and U* as the optimal solution of problem (OH . We 
prove the theorem by contradiction. Assume that U a is not optimal. Then there must exist an e > 0 such that 
U a + 2e < U*. Then there must be a solution x. so that 


U a + 2e < U(xk). 

Let to be sufficiently large so that for any t > to we have 

U(x®) <U a + e. 


Combining ( l33l > and ( OH . we have Uix'^’) + e < U(x.®). 

Let ac be a positive number that satisfies ac < inf |||<9W(xW)|| j, 


for all t. It follows that 


> 

> 

> 

> 


xjj. t+1) - x fe || 2 = ||xj^ - 6®dU® - Xfell 2 

xW-x.ll 2 + (6®) 2 \\dU®\\ 2 - 2 «jW(aw(‘)) ff ( X W -x fe ) 

- (S^) 2 \\du ^\\ 2 - 2<yW(w( x W) - W(x fc )) 




(t) - 112 

Xj. — Xfc 11 


(5 {t) ) 2 ac 2 + 26® e 




- x.ll 2 + 25® e >■■■ 


r (to) 


- Xfc 


2 e^S®. 

j=t 0 


Note that the first inequality is due to the property of subgradient. So we finally have ||x 


(t+i) 


(33) 


(34) 


(35) 


-x.|| 2 >||xW ) - 


Xfe 


+ 2which cannot hold for sufficiently large t. Thus Algorithm [3] optimally solves problem (OH . 


Theorem 3. The optimal solution to problem (ED is also feasible and optima! to problem ( 1291) . 

Proof: From problem (OH to ( [29b . we relax the variables from binary to real and introduce an equality constraint. 
The equality constraint does not change the problem. So the optimal value to problem ( [29b provides an upper bound 
to that of problem (OH . However, it can be observed from Algorithm [3] that the solutions to problem ( [29b are integers 
rather than fractions. So the solutions are also feasible to problem (OH . Since the solutions to problem (OH cannot 
result in a higher optimal value than the solutions to problem ( l29b . the solutions to problem ( [29b are exactly the 
solutions to problem ( OH as well. Henceforth, even though we transform problem ( OH to problem ( l29b . the optimal 
solution is not affected by the transformation. ■ 

To sum up, the optimal solution to problem ( [29b can be solved with Algorithm [3 For comparison purpose, we 
also propose two greedy algorithms as benchmarks, which are presented in Algorithms [4] and 0 The main idea of 
the greedy algorithms is to first identify the most desirable user-BS pair, and then to allocate all the resource to 
that user. This is repeated until convergence is reached. 





Algorithm 3: Two Layer Dual Decomposition Algorithm for Optimization Problem (I3U) 

1 t 

= 0, A (1) = 0 ; 

2 while not converged do 

3 

1 1- t+ 1 ; 

4 

for k = 1, • • • , K do 

5 


for j = 1, • • • , J do 

6 


| Compute Ckj as in 1131 ; 

7 


end 

8 


Find j* = argmaxj |log(c fcj - A®) j ; 

9 


Let = 0 for j ^ j* ; 

10 


if log(cfc 3 . - A^) > 0 then 

11 


1 x k* = 1 i 

12 


else 

13 


= 0 ; 

14 


end 

15 

end 

16 

for j = 1, • • • , J do 

17 


Each BS chooses a step size 0 ' 1 ‘ and computes = min (Lj,eVj 11 } and 

18 

end 

19 end 



Algorithm 4: Greedy Algorithm 4 for Joint Resource Allocation and User Association 

1 Initialize K, = {1, 2, • • • , A'} and J = {1, 2, • • • , J} Xk. to be an all-zero matrix ; 

2 for k = 1 to K do 

3 for j = 1 to J do 

4 | Compute Ckj as in CCO ; 

5 end 

6 end 

7 while maxfcj logics) > 0 do 

8 Find ( k*,j*) = argmax fcj . log(c fcj ) ; 

9 x k * = 1 ; 

3 

10 K. = K.\k* ; 

11 J = J\r ; 

12 end 













Algorithm 5: Greedy Algorithm 5 for Joint Resource Allocation and User Association 


1 Initialize K. = {1, 2, • • • , A'}, J = {1, 2, • • ■ , J} and Xkj to be an all-zero matrix ; 

2 for k = 1 to K do 

3 for j — 1 to J do 

4 | Compute Ckj as in O ; 

5 end 

6 end 

7 for j = 1 to J do 

8 if max fc log(cfe !: ;) > 0 then 

9 Find ( k*,j ) = argmax^c^j ; 

10 Xk * = 1 ; 

j 

11 TC = K.\k* ; 


12 end 

13 end 


IV. Distributed User Association 


In the previous section, we assume a central controller that has global information and assigns users to the BS’s. 
In this section, we consider distributed user association. We still assume that the BS’s have all the CSI via uplink 
training. We further assume that all the BS’s, including the massive MIMO microcell BS and the small cell BS’s, 
belong to the same service provider. Each user makes its own decision based on the broadcast and local information. 
Throughout this section, we do not allow fractional connection. We omit constraint ([2} in the problem formulation, 
which is, however, enforced when solving the problem. 

We model the behavior and interactions among the service provider and users using repeated game theory. The 
first key problem is to determine whether the game will converge. The second key problem is to analyze whether 
both sides are satisfactory about the outcome of the game, i.e., existence of the Nash Equilibrium. 


A. Sen’ice Provider Sets the Price 

The players of the repeated game include the service provider and the users. During each round of the game, the 
service provider determines the price of the connection service. The users decide whether or not to connect, and if 
to connect, to which BS. The strategy of the service provide is to set the price pkj of each BS j for each user k, 
while the strategy of each user k is to set Xk, to either 0 or 1 for j £ J. 

The utility of the service provider is defined as Ub = StLi X//_i x k } Pk :i ■ Since each BS is constrained by its 
maximum load capacity Lj, the service provider aims to solve the following problem. 


K j 



(36) 



k 


The utility of each user is the data rate achieved minus its payment. So each user aims to solve the following 








problem. 


max U k = max < oj k log 



s.t. y jXk . < i, 

3 


E 


x kjPj k i 


J=1 


(37) 


where the logarithmic function represents the satisfaction level of a user k towards its achievable rate, and ui k is 
a weight used to tradeoff rate satisfaction and monetary payment. We assume that the weight u> k of each user is 
drawn from a finite set W with |W| elements. This assumption is true in real-world practice. For instance, $30 
for a wireless service with 60 Mbps data rate is considered to be cheap; $45 is considered to be reasonable; $60 
would be acceptable; $80 would be expensive for most people; $100 would be too expensive; and $150 or above 
would not be an option for most people. So the weight of the users has generally finite choices of values based on 
common sense, and is typically in a range = (0, Wm), where Wm is the maximum possible value for uj k . 

The repeated game is played as follows. Initially, the service provider sets a price for each BS for each user 
and broadcasts the prices to the users. Knowing the prices, the users will feedback the service provider of their 
choices based on their own calculations. Then the service provider updates the prices and broadcasts them to the 
users. Users again inform the service provider of their choices, and so forth. The process is repeated until both the 
service provider and users are all satisfied with the price. 

Given the players, their strategies and utilities, we have the following definition for the NE of the user association 
game. 


Definition 2. A strategy set ^P kj ■ x kj j> f or a ^ j, is an NE of the repeated game ifUB{p * k ., x\.) > UB{p kj , x * k .), 

for all p kj and U k [p* k . ,x kj )> U k {p * kj , x kj ), for all k, x kj . 


Due to the constraint that each user can only connect to one BS, uj k log Xk i Ck i) = l x kj oj k log(cfc 3 -). 

Therefore the objective function of problem (l37l > becomes 


U k = max 


J J 

E Xk ^ k l 0 S(%) - E Xk i p i *’ 

j =i i=i 



(38) 


For the reformulated problem (l38l >. the constraint ■ x kj < 1 indicates that a user may choose not to connect to 
any of the BS’s. On the other hand, if we restrict ^ x kj = 1, then even if the service provider sets the prices to 
infinity, each user will still connect to a BS, which is clearly unreasonable. 

Given the utility function ( |38| > and the constraint in (l37l ). the optimal solution for each user can be derived as 


j* = argmax [u k log (c kj ) - p jk ] (39) 

{ 1 , if j = j* and ut k log(c fe .) > p r 

(40) 

0, otherwise. 

Such users’ decision can be interpreted this way. A user will choose the best connection based on its own 
evaluation. If its evaluation of the connection is greater than or equal to the price, it will connect to this BS. 
Otherwise, the user will not connect to the BS. So we readily have the following result. 


Lemma 5. The highest profit the service provider can obtain from a user k towards BS j, is the user’s evaluation. 

The service provider aims to solve problem (l36l> by tuning variables p kj , k = 1,2 ,K, j = 1,2, ••• , J. 

However, the constraint x kj < Lj is implicitly coupled with all the p k f s, since according to the user’s choice. 


j* = argmax Jg j [oj k log(cfc 3 ) — p ]k \. The service provider problem is actually with the following form. 

K j 

max Ub Xk.pki (41) 

{pkA ,,, 

3 k— 1 3—1 

s.t. x kj( Pk .) — Lji J = 1) 2, • • • , J. 

k 3 

Since problem (|4TT > has coupling constraints, one may try to introduce Lagrange multipliers to the constraint and 
solve the resulting problem using dual decomposition. However, since p kj is implicitly contained in the constraint, 
the gradient and subgradient are difficult to find. Next, we propose Algorithm [6] for the service provider, and then 
prove that the algorithm achieves optimal utility for the service provider and the users. 

Theorem 4. If the service provider adopts Algorithm \6\ the game converges and the NE can be achieved. 

Proof: We first notice that the service provider has priority over the users. The users always make decisions 
based upon the service provider’s price setting. Basically, the service provider controls when the repeated game 
terminates. 

In Algorithmic the service provider tests out the weight of each user using binary search with 0(log 2 (|W|)) 
steps. Once the service provider obtains u>k, k = 1,2, • ■ • , K, it then estimates the users’ price evaluation matrix 
V as follows. 

Vkj = CkjOJk, (42) 

where v kj is the entry of matrix V at row j and column k. Following Lemma 0 the service provider can obtain 
its optimal price strategy by first selecting users for each BS and solving the following problem. 

K J 

max EE XkjVkj (43) 

Xk 3 fc=1 j=1 

s.t. 'y ' Xkj < - Lj , j = 1,2, ■ • ■ , J 

k 

E Xk 3 - 1 ’ k=1 ’ 2 ’--- ,K 

j 

Constraints <[2}. <(42}■ 


The optimal solution x* k . to the above problem can be solved in a similar way as solving problem Pl-2. Then 
the optimal prices for the service provider can be obtained as follows. 


Pk n = 


Vhj, ifx* k .=l; 
v k . + e, otherwise, 


(44) 


where e is an arbitrary positive number. 

Therefore, by adopting Algorithm [6] the optimal utility (highest) can be reached for the service provider. 
Meanwhile, we could see that all the users’ utility must be 0 due to the optimal price setting (i.e., each user’s 
rate satisfaction matches its monetary payment). That means, all the users achieve the optimal utility given the 
price setting as well. Therefore, the game converges to the NE. ■ 

Note that it is possible that the optimal utility of the service provider will be lower than the maximum utility 
during the game, because the load capacity constraint may be violated due to the distributed operation. 


B. A User Bidding based Approach 

We next consider a bidding approach to the problem. Before service starts, users bid to the service provider 
according to their predicted satisfaction towards each BS. And service provider determines whether or not to accept 
a user’s bid and feedback the decisions to users. Then the users make another round of bids according to its predicted 


Algorithm 6: Algorithm for Service Provider 


1 Initialize Umax> u min > t — 0 ; 

2 for k = 1 to K do 

3 for j — 1 to J do 

4 | Compute Ckj as in O ; 

5 end 

6 end 

7 for k = 1 to K do 

8 u k (t ) = ujmax ; 

9 u l k (t) = ujmin ; 

10 end 

n while not converged do 
for k = 1 to K do 

Wk(t) = ^(u k (t) + ; 

for j = 1 to J do 

| Pkjit) = max {wfc(t) logics^ ),0} ; 

end 

end 

t tr- t + 1 ; 

for k = 1 to K do 
if |-Ffc| > 1 then 

1); 

4(f)=4(f-l) ; 

else if | F k \ = 1 then 

ut(t)=ut{t- 1) ; 
w[(t) = wfc(t) ; 

else 

Wfe(t) = Uk(t ) ; 

wUt) =u[{t- 1) ; 

end 


12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 end 

31 end 

32 for /c = 1 to K do 

33 for j = 1 to J do 

34 | Calculate v^ j as in ((42} using Wk ; 

35 end 

36 end 

37 Solve {43]) and find optimal price as in (44[ ; 














Algorithm 7: Algorithm for the Service Provider with the Bidding Approach 

l while not converged do 

2 

for j = 1 to J do 

3 


if BS j is bidden by < Lj users then 

4 


| Keep all the users in BS j’s waiting list ; 

5 


else 

6 


| Keep the top Lj users with the highest bids and reject the other users ; 

7 


end 

8 

end 

9 end 



satisfaction and the service provider’s decision history. The service provider again decides whether or not to accept 
a user’s bid and feedback the decision, and so forth. 

Assume date-intensive users that strive for as high data rate as possible. Each user solves the following problem. 


max Lik = max 

{p kj } 


Yx kj u k log{c ki ),0 
j =l 


S.t. 'Y x kj (p kj ) < 1- 
3 


(45) 


On the other hand, the service provide aims to maximize its utility, i.e., the total payment made by all the users. 


K j 

max Ub = EE Xkj Pkj (46) 

) , , ■ 
o k =1 j=l 

s.t. 'y ' Xkj < Lj, j = 1,2, • ■ • , J. 

h 

Note that the decision variables in these two problems are different from those in problems (l36l > and < 137b . respectively. 

We assume the general case that K > Lj (i.e., not all the users can be served). In order to achieve the 

greatest level of satisfaction, each user makes the highest possible payment. So the optimal solution for each user 
is 


Pkj = max 


I 


x fci w fc log(cfe 3 .),0 


(47) 


The optimal strategy for the service provider is summarized in Algorithm [7] 

During the first stage of the game, each user offers a price to its most desirable BS. Algorithm [7] is used to check 
if each BS j receives more than Lj bids. The service provider only puts Lj top users on BS j’s waiting list based 
on the offered prices; and rejects all other users. If BS j receives no more than Lj bids, all these users will be put 
on BS j’s waiting list. 

At the second stage, if a user is in a BS’s waiting list, it will keep on bidding the same BS with the same price to 
guarantee the highest utility. However, if a user gets rejected in the previous round, as being selfish, it will exclude 
the BS’s that have rejected it and offers a price to its most desirable BS among the remaining ones. For the service 
provider, it adopts the same strategy. If the number of bids received for a BS outnumbers the load capacity of that 
BS, the service provider only keeps the Lj most desirable users on the waiting list and rejects the others. It keeps 
all users on the waiting list if the number of offers received is less than a BS’s load capacity. This two stages repeat 
until convergence is achieved. 


Lemma 6. The sequence of bids made by a user is non-increasing in the user’s preference list. 







Proof: Before a user makes an offer, it computes the satisfaction of all the BS’s to obtain a preference list. 
Since a user aims to maximize its utility, it first proposes to the BS with the highest satisfaction. If it is rejected by 
the BS, it will propose to the BS with the second highest satisfaction, and so forth. Note that even if a user may 
be on the waiting list of a BS, it may be removed from that waiting list at a later stage. If that happens, this user 
will start bidding to other BS. A user will repeat this procedure until it is finally in a BS’s serving list or rejected 
by all BS’s. This concludes the proof. ■ 

Lemma 7. The sequence of bids a BS put on the waiting list is non-decreasing in its preference list. 

Proof: Given the fact any BS has a finite load capacity and I\ > Lj, all the base station will have at 

least one user bidding to it at some stage of the game. Since a BS aims to maximize its utility, it puts all the users 
who make an offer on the waiting list. On the condition that there are too many users, it will reject the users who 
it will never served. In the next round of game, the BS will often have more or at least the same amount of bids 
compared to its current waiting list. This means that the BS has more choices. The BS again only keeps the most 
profitable ones and reject or remove the others from the waiting list. So the sequence of bids a base station put on 
the list is non-decreasing in its preference list. ■ 

Theorem 5. The repeated bidding game converges. 

Proof: Based on Lemmas [6] and [7] we prove this theorem by contradiction. Suppose that this repeated game 
does converge. Then there must be a stage of the game that (i) there is a user k and BS j pair so that user k is 
connected to another BS f or is not connected to any BS; (ii) user k prefers BS j to BS j’ or prefers to be not 
connected; and fiii) BS j prefers user k to a user k! who is on its serving list. 

Consider the case where user k is served by BS f. Since the sequence of bids made by a BS is non-decreasing, 
it must be the case that user k has never bidden to BS j during the game. Otherwise, if user k has bidden to BS 
j, BS j would not have ended up with choosing k! over k. In this case, user k would never have bidden to BS j' 
either, since user k prefers j to j' and the bids (see Lemma |5). However, user k is now served by BS f, user k 
must have bidden to BS j 1 , which contradicts that user k would never have bidden to BS j'. 

The same reasoning holds for the case when user k is not connected to any BS. If BS j prefers k to k' on the 
serving list, BS j would never reject user k while keeping user k!. 

Therefore, the game converges when every user is either on a waiting list or has been rejected by every BS, and 
the game will converge. ■ 

From the proof, we can actually see that the game terminates when the least popular BS becomes fully loaded. 

Theorem 6. The outcome of the repeated bidding game is optimal for both the users and sendee provider. 

Proof: Suppose that the outcome of the game is not optimal for a user k, who is connected to BS j. Then 
there must be another BS j', which has higher ranking than BS j in the preference list of user k and has a serving 

list of users j /[ ,f l . • • • ,j' L l |. Since BS f serves these users, it means that BS f prefers them to user k and BS 

j' is at the top of the preference lists of these users. If at some stage, user k is in the waiting list of BS j (or it is 
inserted by force), the game must have not terminated. 

Since user k is in the waiting list, then one of the final users j[. j' 2 , ■ ■ ■ ,j' L / must be currently off the list, say 
user j' L ( . Then user j' L ; will immediately bid for BS j', since BS j' is at the top of its preference list among all 

the remaining BS’s. And BS j' will remove user k from its waiting list, since user k has a lowest ranking in the 
preference list of BS j'. Thus when the repeated game terminates, the outcomes are optimal for each user. It is 
obvious that the outcome is also optimal for the service provider as well. ■ 

From Theorems 0 and [6j we conclude that the game converges to the NE when the game terminates. 


V. Simulation Study 

We validate the proposed user association schemes with simulations. Throughout the simulations, we assume 
f k = 1/(1 + (%p ) 3 ' 5 ) for the path loss between a user and the massive MIMO BS, and f k = 1/(1 + (ppp) 4 ) for 
the path loss between a user and a small cell BS tm We assume that the power of small scale fading follows a 
uniform distribution from [0.8,1], We fix the location of the massive MIMO BS at the center of the cell. The other 
BS’s are randomly placed across in the cell. Users are randomly placed in the area. The other parameter settings 
are listed in Table [j] The error bars in the plots are 95% confidence intervals. 


TABLE I 

System Configuration 


Parameter 

Value 

Parameter 

Value 

Mmassive 

100 

M 

4 

Lmassive 

10 

L 

4 

p 

- 1 masszve 

40 dBm 

P 

40 dBm 

Area 

1000 X 1000 m 2 

J 

ii 


TABLE II 

Rate Maximization of Centrafized Controe 


K 

50 

100 

150 

200 

250 

Optimal Rate Maximation 

382.9 

483.4 

543.1 

572.3 

594.0 

Greedy Algorithm 1 

363.5 

480.0 

540.4 

571.1 

592.6 

Greedy Algorithm 2 

191.0 

279.6 

340.6 

371.5 

392.5 


TABLE III 




Log Rate Utility of Centralized Control 


K 

50 

100 

150 

200 

250 

Optimal Log Rate Max. 

128.7 

155.5 

167.7 

172.7 

176.0 

Greedy Algorithm 1 

115.0 

153.3 

166.9 

172.4 

175.7 

Greedy Algorithm 2 

67.5 

97.7 

122.0 

133.5 

140.0 


Table [II] presents a comparison of rate maximization with the optimal solution and the two proposed greedy 
algorithms. Tabel [HI] shows a comparison of rate maximization considering proportional fairness with the optimal 
solution and the two proposed greedy algorithms. We can see from both tables that the optimal solution achieves 
the highest network utility. We also notice that as the number of users increases, the gaps between the optimal 
utility and the greedy solutions become more and more narrower. This is because that as there are more users, the 
user diversity effect becomes stronger. So the greedy algorithms and the optimal user association algorithm tend to 
produce similar solutions. 

Throughout this paper, the constraint for each user is Ylj =1 x k, < 1- It should provide upper bounds for the 
problem with the constraint J2j=i x fc. 7 = 1- A comparison of these two different constraints is presented in Fig. Q] 
For a fair comparison, we have exactly the same number of active users as the load capacity for all the BS’s. For 
instance, when the system load capacity is 250, we have J = 51 BS’s and K = 250. We can see that the inequality 
constraint problem indeed upper bounds the equality constraint problem. This is because the inequality constraint 
problem could eliminate the users whose rate is too low with a negative utility. 

Table [TV] presents a comparison of the optimal joint resource allocation and user association algorithm and the 
two proposed greedy algorithms. We find that the optimal scheme achieves the highest utility. Moreover, the gap 
between the optimal scheme and the greedy schemes is quite large. We also consider the equality constraint problem 
as a benchmark for the comparison. For a fair comparison, we set the sum capacity of this system equal to the 
number of users. So there are totally K = 50 active users in the system. The optimal solution of problem P3-2 
achieves a network utility of —59.8462, while the optimal solution of problem d22l > has a network utility of 29.5433. 
We also found that if we connect every user, some edge users will be harmful for the network utility. 

Fig. [2] shows the utility of the service provider and all users when the service provider sets the price (as in 
Section II V- Al l. It can be seen that the repeated game converges after 8 rounds. Furthermore, the utility of all users 
is monotonically decreasing. That is because once a user’s evaluation is known to the service provider, the service 
provider will set prices for the highest profit, which results in 0 utility for that user. As discussed, the utility for 


























Fig. 1. A comparison of the centralized algorithms with logarithmic rate utility under inequality and equality constraints. 


TABLE IV 

Joint Resource Allocation and User Association 


K 

50 

100 

150 

200 

250 

Optimal Joint Resource Allocation 

42.3 

50.8 

55.8 

56.1 

62.0 

and User Association 






Greedy Algorithm 4 

35.3 

37.9 

39.3 

39.9 

40.3 

Greedy Algorithm 5 

35.1 

37.9 

39.3 

39.8 

40.2 



Fig. 2. Convergence of the repeated game when the service provider sets the price and K = 100. 


the service provider is not monotonically increasing, since during the game, the load capacity constraint may be 
violated. Fig. [3 plots the utilities of the service provider and users versus the number of users. We can see that 
as the number of user increases, utility of the service provider also increases. This is mainly due to the effect of 
multi-user diversity. We can also observe that the game terminates after about 8 rounds no matter how many users 
are active. 

Fig-ffldepicts the process of the game when users bid for BS’s (as in Section lTV-Bb . Here we deploy J = 41 BS’s. 




























Fig. 3. Utility of the service provider, utility of the users, and the number of rounds for convergence for systems with various numbers of 
users. 



Fig. 4. Convergence of the repeated game with respect to BS load when users bid. 


The massive MIMO BS has M = 400 antennas. There are I\ = 350 users. The left-hand-side y-axis represents 
the load of the 41 BS’s. The right-hand-side y-axis represents the utility of the service provider. We find the game 
converges in about 10 rounds, and the utility of the service provider is monotonically increasing as the game 
continues. 

To encourage offloading from the macro BS, we consider rate bias for the BS’s in this experiment. Specifically, 
we multiple the rate of the massive MIMO BS with a factor of 0.5 to encourage connection to the PBS’s. Fig. 0 
shows the result when configuration shown in Table |l] It can be observed that the utility with rate bias is higher 
than the utility without considering rate bias. This result demonstrates the efficacy of rate bias and offloading. It 
can be seen that both games terminate in less than 8 rounds. 


VI. Conclusions 

In this paper, we investigated the user association problem in a massive MIMO HetNet from the centralized and 
distributed perspectives. Particularly, by leveraging totally unimodularity we developed optimal algorithms for rate 
maximization and rate maximization with proportional fairness problems. We also developed optimal algorithms 
to the joint resource allocation and user association problem with primal decomposition and dual decomposition. 
Modeling the behavior and interaction of the service provider and users with repeated games, we developed effective 
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Fig. 5. Utility of the service provider with or without rate bias, and convergence of the games under different numbers of users. 


distributed algorithms with proven convergence to the NE. Simulation results verify the efficacy of the proposed 
schemes. 
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